Lead optimization, pharmacophore development and scaffold design of protein kinase CK2 inhibitors as potential COVID-19 therapeutics.

Therapeutic agents being designed against COVID-19 have targeted either the virus directly or the host cellular machinery. A particularly attractive host target is the ubiquitous and constitutively active serine-threonine kinase, Protein kinase CK2 (CK2). CK2 enhances viral protein synthesis by inhibiting the sequestration of host translational machinery as stress granules and assists in viral egression via association with the N-protein at filopodial protrusions of the infected cell. CK2 inhibitors such as Silmitasertib have been proposed as possible therapeutic candidates in COVID-19 infections. The present study aims to optimize Silmitasertib, develop pharmacophore models and design unique scaffolds to modulate CK2. The lead optimization phase involved the generation of compounds structurally similar to Silmitasertib via bioisostere replacement followed by a multi-stage docking approach to identify drug-like candidates. Molecular dynamics (MD) simulations were performed for two promising candidates (ZINC-43206125 and PC-57664175) to estimate their binding stability and interaction. Top scoring candidates from the lead optimization phase were utilized to build ligand-based pharmacophore models. These models were then merged with structure-based pharmacophores (e-pharmacophores) to build a hybrid hypothesis. This hybrid hypothesis was validated against a decoy set and used to screen a diverse kinase inhibitors library to identify favored chemical features in the retrieved actives. These chemical features include; an anion, an aromatic ring and an H-bond acceptor. Based on the knowledge of these features; de-novo scaffold design was carried out which identified phenindiones, carboxylated steroids, macrocycles and peptides as novel scaffolds with the potential to modulate CK2.Communicated by Ramaswamy H. Sarma.

Insights into the structure and dynamics of SARS-CoV-2 spike glycoprotein double mutant L452R-E484Q.

The Receptor Binding Domain (RBD) of SARS-CoV-2, located on the S1 subunit, plays a vital role in the virus binding and its entry into the host cell through angiotensin-converting enzyme 2 (ACE2) receptor. Therefore, understanding the dynamic effects of mutants on the SARS-CoV-2 RBD is essential for discovering drugs to inhibit the virus binding and disrupt its entry into the host cells. A recent study reported a double mutant of SARS-CoV-2, L452R-E484Q, located in the RBD region. Thus, this study employed various computational algorithms and methods to understand the structural impact of both individual variants L452R, E484Q, and the double mutant L452R-E484Q on the native RBD of spike glycoprotein. The effects of the mutations on native RBD structure were predicted by in silico algorithms, which predicted changes in the protein structure and function upon the mutations. Subsequently, molecular dynamics (MD) simulations were employed to understand the conformational stability and functional changes on the RBD upon the mutations. The comparative results of MD simulation parameters displayed that the double mutant induces significant conformational changes in the spike glycoprotein RBD, which may alter its biological functions. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-022-03151-0.

Anti-Fungal Drug Anidulafungin Inhibits SARS-CoV-2 Spike-Induced Syncytia Formation by Targeting ACE2-Spike Protein Interaction.

Drug repositioning continues to be the most effective, practicable possibility to treat COVID-19 patients. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus enters target cells by binding to the ACE2 receptor via its spike (S) glycoprotein. We used molecular docking-based virtual screening approaches to categorize potential antagonists, halting ACE2-spike interactions by utilizing 450 FDA-approved chemical compounds. Three drug candidates (i.e., anidulafungin, lopinavir, and indinavir) were selected, which show high binding affinity toward the ACE2 receptor. The conformational stability of selected docked complexes was analyzed through molecular dynamics (MD) simulations. The MD simulation trajectories were assessed and monitored for ACE2 deviation, residue fluctuation, the radius of gyration, solvent accessible surface area, and free energy landscapes. The inhibitory activities of the selected compounds were eventually tested in-vitro using Vero and HEK-ACE2 cells. Interestingly, besides inhibiting SARS-CoV-2 S glycoprotein induced syncytia formation, anidulafungin and lopinavir also blocked S-pseudotyped particle entry into target cells. Altogether, anidulafungin and lopinavir are ranked the most effective among all the tested drugs against ACE2 receptor-S glycoprotein interaction. Based on these findings, we propose that anidulafungin is a novel potential drug targeting ACE2, which warrants further investigation for COVID-19 treatment.

Targeting cathepsins: A potential link between COVID-19 and associated neurological manifestations.

Many studies have shown that the lysosomal cathepsins, especially cathepsins B/L (CTSB/L) are required for SARS-CoV-2 entry into host cells. Lysosomal proteases, cathepsins are indispensable for normal health and are involved in several brain disorders occurring at different development age periods. On the other hand, it has been well known that COVID-19 infection is largely associated with several neurological disorders. Taken together these findings and given the high levels of expression of CTSB/L in the brain, we here proposed a reasonable hypothesis about the involvement of CTSB/L in the neurological manifestations linked to COVID-19. Pharmacological inhibitions of the CTSB/L could be a potential therapeutic target to block the virus entry as well as to mitigate the brain disorders. To this end, we utilized the network-based drug repurposing analyses to identify the possible drugs that can target CTSB/L. This study identifies the molecules like cyclosporine, phenytoin, and paclitaxel as potential drugs with binding ability to the CTSB/L. Further, we have performed molecular docking and all-atom molecular dynamics (MD) simulations to investigate the stability of CTSL-drug complexes. The results showed strong and stable binding of drugs with CTSL.

Structural stability predictions and molecular dynamics simulations of RBD and HR1 mutations associated with SARS-CoV-2 spike glycoprotein.

The COVID-19 pandemic is caused by human transmission and infection of Severe Acute Respiratory Syndrome Corona Virus-2 (SARS-CoV-2). There is no trusted drug against the virus; hence, efforts are on discovering novel inhibitors for the virus. The entry of a SARS-CoV-2 virus particle into a host cell is initiated by its spike glycoprotein and host Angiotensin-Converting Enzyme 2 (ACE2) receptor interaction. Spike glycoprotein domains, namely, the Receptor Binding Domain (RBD) and Heptad Repeat (HR) domains, are essential for this activity. We have studied the impact of mutations such as A348T, N354D, D364Y, G476S, V483A, S494D in the RBD (319-591), and S939F, S940T, T941A, S943P (912-984) in the HR1 domains of spike glycoprotein. Summarily, we utilized the computational screening algorithms to rank the deleterious, damaging and disease-associated spike glycoprotein mutations. Subsequently, to understand the changes in conformation, flexibility and function of the spike glycoprotein mutants, Molecular Dynamics (MD) simulations were performed. The computational predictions and analysis of the MD trajectories suggest that the RBD and HR1 mutations induce significant phenotypic effects on the pre-binding spike glycoprotein structure, which are presumably consequential to its binding to the receptor and provides lead to design inhibitors against the binding.Communicated by Ramaswamy H. Sarma.

Simultaneous Inhibition of SARS-CoV-2 Entry Pathways by Cyclosporine.

The COVID-19 pandemic caused by SARS-CoV-2 represents a global public health emergency. The entry of SARS-CoV-2 into host cells requires the activation of its spike protein by host cell proteases. The serine protease, TMPRSS2, and cysteine proteases, Cathepsins B/L, activate spike protein and enable SARS-CoV-2 entry to the host cell through two completely different and independent pathways. Therefore, inhibiting either TMPRSS2 or cathepsin B/L may not sufficiently block the virus entry. We here hypothesized that simultaneous targeting of both the entry pathways would be more efficient to block the virus entry rather than targeting the entry pathways individually. To this end, we utilized the network-based drug repurposing analyses to identify the possible common drugs that can target both the entry pathways. This study, for the first time, reports the molecules like cyclosporine, calcitriol, and estradiol as candidate drugs with the binding ability to the host proteases, TMPRSS2, and cathepsin B/L. Next, we analyzed drug-gene and gene-gene interaction networks using 332 human targets of SARS-CoV-2 proteins. The network results indicate that, out of 332 human proteins, cyclosporine interacts with 216 (65%) proteins. Furthermore, we performed molecular docking and all-atom molecular dynamics (MD) simulations to explore the binding of drug with TMPRSS2 and cathepsin L. The molecular docking and MD simulation results showed strong and stable binding of cyclosporine A (CsA) with TMPRSS2 and CTSL genes. The above results indicate cyclosporine as a potential drug molecule, as apart from interacting with SARS-CoV-2 entry receptors, it also interacts with most of SARS-CoV-2 target host genes; thus it could potentially interfere with functions of SARS-CoV-2 proteins in human cells. We here also suggest that these antiviral drugs alone or in combination can simultaneously target both the entry pathways and thus can be considered as a potential treatment option for COVID-19.

Screening Malaria-box compounds to identify potential inhibitors against SARS-CoV-2 Mpro, using molecular docking and dynamics simulation studies.

Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) Main protease (Mpro) is one of the vital drug targets amongst all the coronaviruses, as the protein is indispensable for virus replication. The study aimed to identify promising lead molecules against Mpro enzyme through virtual screening of Malaria Venture (MMV) Malaria Box (MB) comprising of 400 experimentally proven compounds. The binding affinities were studied using virtual screening based molecular docking, which revealed five molecules having the highest affinity scores compared to the reference molecules. Utilizing the established 3D structure of Mpro the binding affinity conformations of the docked complexes were studied by Molecular Dynamics (MD) simulations. The MD simulation trajectories were analysed to monitor protein deviation, relative fluctuation, atomic gyration, compactness covariance, residue-residue map and free energy landscapes. Based on the present study outcome, we propose three Malaria_box (MB) compounds, namely, MB_241, MB_250 and MB_266 to be the best lead compounds against Mpro activity. The compounds may be evaluated for their inhibitory activities using experimental techniques.

Targeting SARS-CoV-2 nucleocapsid oligomerization: Insights from molecular docking and molecular dynamics simulations.

The outbreak of COVID-19 caused by SARS-CoV-2 virus continually led to infect a large population worldwide. Currently, there is no specific viral protein-targeted therapeutics. The Nucleocapsid (N) protein of the SARS-CoV-2 virus is necessary for viral RNA replication and transcription. The C-terminal domain of N protein (CTD) involves in the self-assembly of N protein into a filament that is packaged into new virions. In this study, the CTD (PDB ID: 6WJI) was targeted for the identification of possible inhibitors of oligomerization of N protein. Herein, multiple computational approaches were employed to explore the potential mechanisms of binding and inhibitor activity of five antiviral drugs toward CTD. The five anti-N drugs studied in this work are 4E1RCat, Silmitasertib, TMCB, Sapanisertib, and Rapamycin. Among the five drugs, 4E1RCat displayed highest binding affinity (-10.95 kcal/mol), followed by rapamycin (-8.91 kcal/mol), silmitasertib (-7.89 kcal/mol), TMCB (-7.05 kcal/mol), and sapanisertib (-6.14 kcal/mol). Subsequently, stability and dynamics of the protein-drug complex were examined with molecular dynamics (MD) simulations. Overall, drug binding increases the stability of the complex with maximum stability observed in the case of 4E1RCat. The CTD-drug complex systems behave differently in terms of the free energy landscape and showed differences in population distribution. Overall, the MD simulation parameters like RMSD, RMSF, Rg, hydrogen bonds analysis, PCA, FEL, and DCCM analysis indicated that 4E1RCat and TMCB complexes were more stable as compared to silmitasertib and sapanisertib and thus could act as effective drug compounds against CTD.Communicated by Ramaswamy H. Sarma.

Insights into the structural and dynamical changes of spike glycoprotein mutations associated with SARS-CoV-2 host receptor binding

Novel Coronavirus or SARS-CoV-2 has received worldwide attention due to the COVID-19 pandemic, which originated in Wuhan, China leading to thousands of deaths to date. The SARS-CoV-2 Spike glycoprotein protein is one of the main focus of COVID-19 related research as it is a structural protein that facilitates its attachment, entry, and infection to the host cells. We have focused our work on mutations in two of the several functional domains in the virus spike glycoprotein, namely, receptor-binding domain (RBD) and heptad repeat 1 (HR1) domain. These domains are majorly responsible for the stability of spike glycoprotein and play a key role in the host cell attachment and infection. In our study, several mutations like R408I, L455Y, F486L, Q493N, Q498Y, N501T of RBD (319-591), and A930V, D936Y of HR1 (912-984) have been studied to examine its role on the spike glycoprotein native structure. Comparisons of MD simulations in the WT and mutants revealed a significant de-stabilization effect of the mutations on RBD and HR1 domains. We have investigated the impact of mapped mutations on the stability of the spike glycoprotein, before binding to the receptor, which may be consequential to its binding properties to the receptor and other ligands. Communicated by Ramaswamy H. Sarma.